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Abstract. We have measured the correlation of rates 
in IceCube with long and short term variations in 
the South Pole atmosphere. The yearly temperature 
variation in the middle stratosphere (30-60 hPa) is 
highly correlated with the high energy muon rate 
observed deep in the ice, and causes a ±10% seasonal 
modulation in the event rate. The counting rates of 
the surface detectors, which are due to secondary 
particles of relatively low energy (muons, electrons 
and photons), have a negative correlation with tem- 
peratures in the lower layers of the stratosphere 
(40-80 hPa), and are modulated at a level of ±5%. 
The region of the atmosphere between pressure levels 
20-120 hPa, where the first cosmic ray interactions 
occur and the produced pions/kaons interact or decay 
to muons, is the Antarctic ozone layer. The anti- 
correlation between surface and deep ice trigger 
rates reflects the properties of pion/kaon decay and 
interaction as the density of the stratospheric ozone 
layer changes. Therefore, IceCube closely probes the 
ozone hole dynamics, and the temporal behavior of 
the stratospheric temperatures. 
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I. Introduction 

The IceCube Neutrino Observatory, located at the 
geographical South Pole (altitude 2835m), has been 
growing incrementally in size since 2005, surrounding 
its predecessor AMANDA] 1]. As of March 2009, Ice- 
Cube consists of 59 strings in the Antarctic ice, and 
59 stations of the IceTop cosmic ray air shower array 
on the surface. Each IceCube string consists of 60 
Digital Optical Modules (DOMs) deployed at depths of 
1450-2450m, and each IceTop station comprises 2 ice 
Cherenkov tanks with 2 DOMs in each tank. 

IceCube records the rate of atmospheric muons with 
E M > 400 GeV. Muon events pass the IceCube Simple 
Majority Trigger (InIce_SMT) if 8 or more DOMs are 
triggered in 5/isec. The IceTop array records air showers 
with primary cosmic ray energies above 300 TeV. Air 
showers pass the IceTop Simple Majority Trigger (Ice- 
Top_SMT) if 6 or more surface DOMs trigger within 
5/isec. 



As the Antarctic atmosphere goes through seasonal 
changes, the characteristics of the cosmic ray interac- 
tions in the atmosphere follow these variations. When 
the primary cosmic ray interacts with atmospheric nu- 
clei, the pions and kaons produced at the early interac- 
tions mainly determine the nature of the air shower. For 
the cosmic ray energies discussed here, these interactions 
happen in the ozone layer (20-120 hPa pressure layer 
at 26 down to 14 km) of the South Pole stratosphere. 
During the austral winter, when the stratosphere is cold 
and dense, the charged mesons are more likely to interact 
and produce secondary low energy particles. During the 
austral summer when the warm atmosphere expands and 
becomes less dense, the mesons more often decay rather 
than interact. 

Figure 1 demonstrates the modulation of rates in 
relation to the temporal changes of the South Pole strato- 
sphere. The scaler rate of a single IceTop DOM on the 
surface is mostly due to low energy secondary particles 
(MeV electrons and gammas, ~1 GeV muons)||2| pro- 
duced throughout the atmosphere, and therefore highly 
modulated by the atmospheric pressure. However, after 
correcting for the barometric effect, the IceTop DOM 
counting rate also reflects the initial pion interactions in 
the middle stratosphere. The high energy muon rate in 
the deep ice, on the other hand, directly traces the decay 
characteristics of high energy pions in higher layers 
of the stratosphere. The IceCube muon rate reaches its 
maximum at the end of January when the stratosphere 
is warmest and most tenuous. Around the same time 
Icetop DOMs measure the lowest rates on the surface 
as the pion interaction probability reaches its minimum. 
The high energy muon rate starts to decline as the 
stratosphere gets colder and denser. The pions will 
interact more often than decay, yielding the maximum 
rate in IceTop and the minimum muon rate in deep ice 
at the end of July. 

II. Atmospheric effects on the IceCube Rates 

The Antarctic atmosphere is closely monitored by the 
NOAA Polar Orbiting Environmental Satellites (POES) 
and by the radiosonde balloon launches of the South 
Pole Meteorology Office. However, stratospheric data 
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Fig. 1. The temporal behavior of the South Pole stratosphere from May 2007 to April 2009 is compared to IceTop DOM counting rate and 
the high energy muon rate in the deep ice. (a) The temperature profiles of the stratosphere at pressure layers from 20 hPa to 100 hPa where 
the first cosmic ray interactions happen, (b) The IceTop DOM counting rate (black -observed, blue -after barometric correction) and the surface 
pressure (orange), (c) The IceCube muon trigger rate and the calculated effective temperature (red). 



is sparse during the winter when the balloons do not 
reach high altitudes, and satellite based soundings fail 
to return reliable data. For such periods NOAA derives 
temperatures from their models. We utilize both the 
ground-based data and satellite measurements/models 
for our analysis. 

A. Barometric effect 

In first order approximation, the simple correlation 
between log of rate change A{lnR} and the surface 
pressure change AP is 



A{lnR} — /3 ■ AP 



(1) 



where f3 is the barometric coefficient. 

As shown by the black line in the Figure lb, the 
observed IceTop DOM counting rate varies by ±10% in 
anti-correlation with surface pressure, and the barometric 
coefficient is determined to be f3 = — 0.42%/hPa. Using 
this value, the pressure corrected scaler rate is plotted 
as the smoother line (blue) in Figure lb. The cosmic 
ray shower rate detected by the IceTop array also varies 
by ±17% in anti-correlation with surface pressure, and 
can be corrected with a f3 value of — 0.77%/hPa. As 
expected [3], the IceCube muon rate shown in Figure 
lc is not correlated with surface pressure. However, 
during exceptional stratospheric temperature changes, 
the second order temperature effect on pressure becomes 
large enough to cause anti-correlation of the high energy 
muon rate with the barometric pressure. During such 



events the effect directly reflects sudden stratospheric 
density changes, specifically in the ozone layer. 

B. Seasonal Temperature Modulation 

Figure 1 clearly demonstrates the seasonal temper- 
ature effect on the rates. The IceTop DOM counting 
rate, after barometric correction, shows ±5% negative 
temperature correlation. On the other hand, the IceCube 
muon rate is positively correlated with ±10% seasonal 
variation. 

From the phenomenological studies (4) 0, it is known 
that correlation between temperature and muon intensity 
can be described by the effective temperature T e ff, 
defined by the weighted average of temperatures from 
the surface to the top of the atmosphere. T e ff approxi- 
mates the atmosphere as an isothermal body, weighting 
each pressure layer according to its relevance to muon 
production in atmosphere [5||6|. 

The variation of muon rate AR^/ < R^ > is related 
to the effective temperature as 



AR„ 



AT, 



< Rn> 



CtT- 



eff 



(2) 



<Teff>' 

where clt is the atmospheric temperature coefficient. 

Using balloon and satellite data for the South Pole 
atmosphere, we calculated the effective temperature as 
the red line in Figure lc. We see that it traces the 
IceCube muon rate remarkably well. The calculated 
temperature coefficient = 0.9 for the IceCube muon 
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TABLE I 

Temperature and correlation coefficients of rates for different 
stratospheric layers of 20-100 hPa and T e ff- 



iceCube muon rate vs. T p (p=40-50hPa} 
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20-30 


214.0 


0.512 


0.953 


-0.194 


-0.834 


30-40 


208.7 


0.550 


0.986 


-0.216 


-0.906 


40-50 


207.3 


0.591 


0.993 


-0.240 


-0.946 


50-60 


206.6 


0.627 


0.985 


-0.261 


-0.968 


60-70 


206.3 


0.656 


0.971 


-0.278 


-0.975 


70-80 


206.3 


0.679 


0.954 


-0.292 


-0.975 


80-100 


206.5 


0.708 


0.927 


-0.310 


-0.971 


<T etf > 


CUT 


7 


CtT 


7 


211.3 


0.901 


0.990 


-0.360 


-0.969 



rate agrees well with the expectations of models as well 
as with other experimental measurements[7||8|. 

In this paper, we also study in detail the relation be- 
tween rates and stratospheric temperatures for different 
pressure layers from 20 hPa to 100 hPa as 

AR P ATP 

= O&r . (3) 

< R > t <Tp> w 
The temperature coefficients for each pressure layer 
Oj, and the correlation coefficient 7 are determined from 
regression analysis. Pressure-corrected IceTop DOM 
counting rate, and IceCube muon rate are sorted in 
bins of ~ 10 days, and deviations AR/ < R > from 
the average values are compared with the deviation of 
temperatures at different depths AT P / < T p >. 

We list the values of ctj, and 7 for IceCube muon 
rate and IceTop DOM counting rate in Table 1. We 
find that the IceCube muon rate correlates best with 
the temperatures of 30-60 hPa pressure layers, while the 
IceTop DOM counting rate shows the best correlation 
with 60-80 hPa layers. In Figure 2. we plot the rate and 
temperature correlation for layers which yield the best 
correlation. 

III. Exceptional Stratospheric Events and 
the Muon rates 

The South Pole atmosphere is unique because of 
the polar vortex. In winter a large-scale counter clock- 
wise flowing cyclone forms over the entire continent 
of Antarctica, isolating the Antarctic atmosphere from 
higher latitudes. Stable heat loss due to radiative cooling 
continues until August without much disruption, and the 
powerful Antarctic vortex persists until the sunrise in 
September. As warm air rushes in, the vortex loses its 
strength, shrinks in size, and sometimes completely dis- 
appears in austral summer. The density profile inside the 
vortex changes abruptly during the sudden stratospheric 
warming events, which eventually may cause the vortex 
collapse. The ozone depleted layer at 14-21 km altitude 
(ozone hole), observed in September/October period, is 
usually replaced with ozone rich layer at 18-30 km soon 
after the vortex breaks up. 
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Fig. 2. Correlation of IceCube muon and IceTop DOM counting rates 
with stratospheric temperatures and T £ jj. (a) IceCube muon rate vs. 
temperature at 40-50 hPa pressure layer, (b) IceCube muon rate vs. 
effective temperature, (c) IceTop DOM counting rate vs. temperature at 
70-80 hPa pressure layer, (d) IceTop DOM counting rate vs. effective 
temperature. 



Apart from the slow seasonal temperature variations, 
IceCube also probes the atmospheric density changes 
due to the polar vortex dynamics and vigorous strato- 
spheric temperature changes on time scales as short as 
days or even hours, which are of great meteorological 
interest. 

An exceptional and so far unique stratospheric event 
has already been observed in muon data taken with 
IceCube's predecessor AMANDA-II. 

A. 2002 Antarctic ozone hole split detected by AMANDA 

In late September 2002 the Antarctic stratosphere 
underwent its first recorded major Sudden Stratospheric 
Warming (SSW), during which the atmospheric temper- 
atures increased by 40 to 60 K in less than a week. 
This unprecedented event caused the polar vortex and 
the ozone hole, normally centered above the South Pole, 
to split into two smaller, separate off-center parts (Figure 
3) 0. 

Figure 4 shows the stratospheric temperatures between 
September and October 2002 along with the AMANDA- 
II muon rate. The muon rate traces temperature varia- 
tions in the atmosphere in great detail, with the strongest 
correlation observed for the 40-50 hPa layer. 

B. South Pole atmosphere 2007-2008 

Unlike in 2002, the stratospheric conditions over 
Antarctica were closer to average in 2007 and 2008. 
In 2007 the polar vortex was off-center from the South 
Pole during most of September and October, resulting in 
greater heat flux into the vortex, which decreased rapidly 
in size. When it moved back over the colder Pole region 
in early November it gained strength and persisted until 
the beginning of December. The 2008 polar vortex was 
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Fig. 3. Ozone concentration over the southern hemisphere on 
September 20th 2002 (left) and September 25th 2002 (right) (TO). 
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Fig. 4. Average temperatures in various atmospheric layers over the 
South Pole (top) and deep ice muon rate recorded with the AMANDA- 
II detector (bottom) during the Antarctic ozone hole split of September- 
October 2002. 



GDAS-CPC Zonal Temperature Time Series 




GDAS-CPC Zonal Temperature Time Series 




Fig. 5. The temperature time series of the Antarctic atmosphere 
produced by NOAAQTj are shown for 2007 and 2008. The pattern 
observed in the deep ice muon rate (black line) is superposed onto 
the plot to display the striking correlation with the stratospheric 
temperature anomalies. 



one of the largest and strongest observed in the last 10 
years over the South Pole. Because of this, the heat flux 
entering the area was delayed by 20 days. The 2008 
vortex broke the record in longevity by persisting well 
into mid-December. 

In Figure 5 we overlay the IceCube muon rate over 
the temperature profiles of the Antarctic atmosphere pro- 
duced by the NOAA Stratospheric analysis team lfTTII . We 
note that the anomalous muon rates (see, for example, 
the sudden increase by 3% on 6 August 2008) observed 
by IceCube are in striking correlation with the middle 
and lower stratospheric temperature anomalies. 

We are establishing automated detection methods for 
anomalous events in the South Pole atmosphere as well 
as a detailed understanding and better modeling of 
cosmic ray interactions during such stratospheric events. 
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